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We study chaos in the Hamiltonian Mean Field model (HMF), a system with many 
degrees of freedom in which A'^ classical rotators are fully coupled. We review the most 
important results on the dynamics and the thermodynamics of the HMF, and in particular 
we focus on the chaotic properties. We study the Lyapunov exponents and the Kolmogorov- 
Sinai entropy, namely their dependence on the number of degrees of freedom and on energy 
density, both for the ferromagnetic and the antiferromagnetic case. 

§1. Introduction 

In systems with a few degrees of freedom the Largest Lyapunov Exponent (LLE), 
which quantifies chaotic motion, is often studied as a function of the control param- 
eter. In many-degrees-of-freedom systems, this can also be done and, moreover, the 
control parameter may acquire a more transparent physical meaning, making refer- 
ence to a thermodynamic quantity. There have been indeed several studiesu^'cP of 
the dependence of the LLE on energy density in Hamiltonian systems with short- 
range interactions (e.g. FPU lattices), which, up to now confirm the conjecture&'S* 
that the LLE reaches a finite, energy dependent, value in the thermodynamic (large- 
volume) limit. Moreover, a scaling limit exists for the full Lyapunov spectrum, which 
implies that the Kolmogorov-Sinai entropy Sks scales with the volume. 

Both the question of the existence of a well defined thermodynamic limit of LLE 
and Sks^ and their dependence on energy (or other control parameters) are open 
for systems with long-range interactions. In this paper we review the most recent 
results on this subject for a Hamiltonian model with many degrees of freedom, named 
Hamiltonian Mean Field (HMF), which describes a, fully- coupled system of N classical 
spins frotators) in the attractive (ferromagnetic) and repulsive (antiferromagnetic) 
casesi^'&. HMF has been recentlv thoroughly investigated both from a theoretical 
and a numerical point of viewQ^'i^'M'S), revealing a deep link between dynamics and 
thermodynamics. Here, we discuss this relation by studying the Lyapunov exponents 
and the Kolmogorov-Sinai entropy, namely their dependence on N and on energy, 
both for the ferromagnetic and the antiferromagnetic case. 
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§2. The model 

The model describes a system of N classical spins (rotators) nij = (cosOi, sinOi). 
Each spin i is characterized by the angle < 0j < 27r and the conjugate momentum 
Pi, and is fully coupled to all the others. The Hamiltonian is: 



H{e,p)=K + V, (2-1) 



where 



TV 2 N 

^ = Ei- .^ = ^ E[i-c«K^^-^.)] (2-2) 

are the kinetic and potential energy. The potential energy V corresponds to the 
interaction of the X — Y model in the infinite-range mean field case, and this is 
the reason why the model has been named Hamiltonian Mean Field (HMF). The 
case e = 1 describes a ferromagnetic behavior, while e = — 1 corresponds to an 
antiferromagnetic interaction. The model has a possible alternative interpretation. 
It can in fact be seen as a system of particles moving on a circle, the position of each 
particle being given by the angle 9i and its momentum by pi. 

The success of the Hamiltonian mean field model is based on the fact that both 
its statistical mechanics and its dynamics can be treated in relatively simple way. 
In fact the thermodynamics of the HMF can be derived exactly for A'^ — oo in the 
canonical ensemble, both for the ferromagnetic and for the antiferromagnetic caseB-* . 
A total magnetization vector can be defined as M = Y^f=i The ferromagnetic 
system has a second-order phase transition from a clustered phase with M = [M| 7^ 
to a disordered phase with M = as a function of energy or temperature. In the 
antiferromagnetic case spins tend to be opposite to each other (interaction among 
the particles is repulsive) and therefore M = (disordered state) for any value of 
the temperature. 

On the other side the dynamics of the system can be investigated for a relatively 
large value of N (we have considered N up to a value of 20000) by solving the 2N 
coupled equations of motion: 

Oi=Pi, Pi = -eMsin{0i-(j)) , i = l,...,N , (2-3) 

where (M, (p) are respectively the modulus and the phase of the total magnetiza- 
tion vector M. Each spin moves in a mean field which is in turn generated self- 
consistently by the all the other spins. In the N ^ 00 limit the dynamics of the 
HMF can be seen as the interaction of a single spin with a mean field, and the equa- 
tions are formally equivalent to those of a perturbed pendulum. Solving equations 



(2-3) corresponds to treating the system in the microcanonical ensemble, because the 
total energy is conserved along each dynamical trajectory (also total momentum is 
a conserved quantity and is typically fixed at zero). 

Hence, HMF has a very remarkable property: it is possible to compare the results 
of the canonical and microcanonical ensemble. Moreover the HMF is a high dimen- 
sional system with long-range forces where one can explore deviations from standard 
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thermodynamics ^S-i)' . HMF exhibits a very rich non-equihbrium dynam- 
ics, and by using the dynamical approach we have studied in detail the problem of 
the relaxation to the canonical equilibrium. The following is a brief review of the 
main results. 

Ensemble inequivalence: because of the long-range nature of the interaction, 
the microcanonical ensemble gives different predictions from the canonical ensemble 
This is true both in the ferromagnetic case, where we fjud the presence of 
quasistationary states different from the canonical equilibriumEF'EP'liy), and in the 
antiferromagnetic case, where a very particular collective phenomenon (bi-cluster 
rmation) appears at low energies in disagreement with the canonical predictions 
I). 

Metastahility: in the ferromagnetic case microcanonical simulations show ±he 
presence of quasi-stationary metastable states with negative specific heatla^'EP'liSl'. 
In fact, if the system is started in far-off-equilibrium initial conditions (for example 
in a "water bag", i.e. putting all the rotators at = and giving them a uni- 
form distribution of velocities with a finite width centered around zero), it does not 
relax directly to the canonical equilibrium. Instead we observe a stabilization into 
metastable states. The temperature of these states are different from the canonical 
one and the velocity distributions are not Gaussianll3. The metastable states are 
called quasistationary states because they have a lifetime which increases linearly 
with the number of particles A^. Thev are expected to become real equilibrium 
solutions in the thermodynamic limit 

Collective phenomena: in the antiferromagnetic bi-cluster of rotators at 

a distance tl in angle is present at very low energy in the microcanonical numerical 
simulations . In fact, below a threshold energy, particles groups spontaneously into 
two big clusters and oscillate maintaning the total magnetization equal to zero. This 
is a pure microcanonical result, stable also for very large and not in agreement 
with the canonical ensemble where a disordered state with all the spins randomly 
oriented is predicted. This collective phenomenon modifies the energy-temnerature 
relation at very small energies and is an effect of the long-range interactionllZP . 

Anomalous diffusion: diffusion and transport of a particle in a medium or in a 
fluid flow are characterized by the average square displacement cr^(t). In general one 
has 

a'^{t) ~ r (2-4) 

with a = 1 for riorrnal diffusion. All the processes with a 7^ 1 are termed anomalous 
diffusion. EJ^'Ej^'EJ'EHP'E^^'EJ) In our model the variance of the spin angle 6 can be 
defined according to the expression: 

al{t) =< {9- <e>f> (2-5) 

where < . > stands for an average over the N spins. Superdiffusion with a = 
1.38 lb 0.05 is observed in the ferromagnetic case injjie energy range 0.5 <U < 0.75, 
i.e. slightly below the critical energy Uc = 0.75I12P. Superdiffusion is due to the 
presence of Levy flights and after a transient regime a change to the slope a = 1 
(normal diffusion) is observedil' . Normal diffusion occurs at a crossover time which 
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we have found to coincide with the relaxation time to canonical equilibrium0 . 
However in the continuum limit diffusion is always anomalous since the crossover 
time and the relaxation one diverge with A^. These irivestigation confirms pioneering 
studies by Kaneko and Konishi for coupled mapsL^^'cH*, and are on the same line of 
investigation of refs.EZP'Es) where the effect of noise or fluctuations on the diffusion 
is studied. 

Generalization of HMF: generalizations of the model have already appeared in 
the literature. In the model introduced by Anteneodo and Tsallis the rotators have 
been attached to the sites of a ID lattice and the interaction of the HMF has been 
modulated by a term dependinff on the lattice distance between two spins, going like 
a decaying power-law r~"EiP'll3'. The 2D case of the HMF has been considered by 
Antoni and TorciniEM). Further generalizations, like ID models with spatial mod- 
ulations or ID models with mixed (attractive-repulsive) interactions, are currently 
under investigation. 

In this paper we focus mainly on the study of the chaotic dynamics in the 
HMF. Our investigations are relevant for the foundation of statistical mechanics 
and also for the study of phase transitions in finite size systems like for example: 
nuclear multifragmentationS-*'® , atomic clustersE!''© and astrophysics©. HMF, 
with the possibility of considering both the ferromagnetic and the antiferromagnetic 
cases, offers two different scenarios. In the following sections we will show that the 
chaotic properties of the system are different in these two cases. 



§3. Thermodynamics: the canonical solution 



The HMF presents the noticeable advantage of possessing an exact solution 
in the canonical ensemble. Therefore microscopic dynamics, and in particular the 
chaotic properties, can be studied in connection with the thermodynamic macro- 
scopic behavior. 

In this section we discuss the canonical solution for both e = +1 (ferromagnet) 
and e = — 1 (antiferromagnet). 

In the ferromagnetic model the potential is attractive and the ground state of the 
system is reached at [/ = where all spins are parallel (all the rotators have the same 
position on the circle). In the antiferromagnetic model the potential is repulsive and 
the ground state, reached aX U = —1/2, consists in a randomly uniform distribution 
of the spins orientations. 

In the high temperature region, both in the ferromagnetic and in the antiferro- 
magnetic model, the rotators are randomly distributed on the circle; each rotator 
moves uniformly around the circle and the modulus of |M| is equal to zero. Ther- 
modynamically, when the potential is attractive the HMF has a second-order phase 
transition with order parameter M, while in the repulsive case the free energy is 
smooth and M = for any value of T. The exact solution of the model in the 
canonical ensemble predicts a caloric curve given by 



2 2^ ' 



(3-1) 
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Fig. 1. We plot the behavior of temperature and magnetization as a function of the energy per 
particle U in the ferromagnetic (a)(b) and in the antiferromagnetic case (c)(d). Points are 
microcanonical numerical simulations for A'' = 100, 1000 while the solid curves are the exact 
canonical predictions, see text. The dashed line indicates the critical point for the ferromagnetic 
case. 



The ferromagnetic case has a critical temperature Tc = 0.5 (Uc = 0.75), while in 

the antiferromagnetic case there is no phase transition and M = in equation ( |3-1| ). 

In fig.l we plot the behavior of magnetization and temperature as a function of the 

energy per particle. The solid curves are the canonical predictions (see eq.(3-l)) while 

the points are the results of the microcanonical numerical simulations for a system 

with N=100 and N=1000 (see next section). The vertical dashed line indicates the 

critical point for the ferromagnetic case. The figure shows that there are small 

deviations due to the finite size of the system. These deviations decrease as N~^^'^. 

In the repulsive case there are also deviations from the canonical equilibrium at very 

small energies due to the formation of the bi-cluster (for a detailed discussion see 
I I. ^ ^ 

ref.F 
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§4. Dynamics: the microcanonical simulations 

We have integrated equations ( |2-3D on the computer by means of a fourth order 
symplectic algorithmo* , with a time step fixed in order to have a good energy 
conservation (the average error is of the order = 10~^, but at very low energy it 
is necessary to have a higher accuracy and simulations were done with = 10~^^). 
For each dynamical trajectory we start the system with a given initial distribution 
and we compute 9i,pi at each time step, and from them the total magnetization 
M and the temperature T (through the relation T = 2 < K > /N). We consider 
systems with various sizes N and we explore a wide range of energies U = E/N. 
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Fig. 2. We show the largest Lyapunov exponent Ai and the Kolmogorov-Sinai entropy per particle 
Sks/N as a function of U — E/N for the ferromagnetic and the antiferromagnetic case. Points 
are numerical simulations for N=100, lines are only to guide the eye. The dashed vertical line 
indicates the critical point. See text. 

To characterize chaos we compute the entire Lvapunov spectrum = 1,2,... , 2N 
by using the standard method of Benettin et al.l^'t^. We focus our attention on 
the Largest Lyapunov Exponent (LLE) Ai and on the Kolmogorov-Sinai entropy 
Sks = J2iLi -^j) computed as the sum of the positive Lyapunov exponentsE3) In fig.2 
we report Ai and Sks/^ vs. U , both for the ferromagnetic and the antiferromagnetic 
case (N=100). The HMF is integrable in the two limits of very low and very high 
energy, and this is indifferently true for e = ±1. The LLE and the Kolmogorov-Sinai 
entropy tend to zero in these two limits. 
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In particular the HMF for e = 1 has aheady been studied in detail and the 
following results have been reported in the literature: 

- for [/ ^ 0, Ai ^ as where the exponent is found to be ^ = 1/2. Essentially 
no dependence on the system size is observed in this regime Q^'EP. 

- for U » Uc, Ai ^ as A^~3 and this behavior can be explained by means of a 
random matrix approximation!^' 0. 

The main difference between the ferromagnetic and the antiferromagnetic hamilto- 
nian appears at intermediate energy. In fact, although both cases are chaotic and we 
have Ai and Sks/^ different from zero, in the ferromagnetic system one observes a 
well defined peak just below the critical energy Uc- Jn some sense, the dynamics feels 
the presence of the phase transition. In fact in ref . B* the increase of the LLE for the 
ferromagnetic case has been related to the increase of kinetic energy fluctuations and 
the specific heat. It has also been shown that the peak, persists as ^ ooi). This 
result is also confirmed by a recent niDre sophisticated theoretical calculation 
using the formalism introduced in ref. . 
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Fig. 3. We show the behavior of Ai vs. A'^ for different energies in the ferromagnetic (a) and in the 
antiferromagnetic case (b). See text for more details 

On the contrary, in the antiferromagnetic case a smoothed shoulder is found 
(instead of a peak) both for Ai and for Sks/^- The difference between the two 
cases is better visible for Sks/^- In a pioneering paper a similar pronouncedDeak 
in LLE was found for second-order phase transitions in nuclear-like systems c3 . A 
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smooth behavior similar to the one in fig. 2(b) was found in other models, when there 
is no phase transition in the canonical ensemble a^'E3^. 

This different behavior can be better shown by studying the LLE as a function 
of the size of the system. This is displayed in figure 3, where we report Ai as a 
function of N for several energies below and above the peak. In the ferromagnetic 
case Ai is constant or even increases for U < Uc, while a decay to zero as N~^^^ 
is evident for U > Uc- In fact, for U > Uc the rotators move independently and 
the power law Aj^.^cx A^~^/'^ can be well explained by a suitable random matrix 
approximation eP' £3). On the contrary for the antiferromagnetic case fig. 3(b), the 
LLE appears to vanish as — > 00 for all values of U (apart from the very low 
energy region, where the bi-cluster forms). In particular, we find the same power 
law Ai oc A^~^/^ as for the ferromagnetic case in the overcritical region. In fact, 
the random matrix approximation applies also when the potential is repulsive. The 
chaotic behavior is only an artifact of the finite-size fluctuations which disappear 
for N 00. No chaos exist in the thermodynamic limit. This behavior strongly 
contrasts with what happens for FPU lattices (short-range interactions), where the 
LLE reaches a finite value, and therefore chaos persists, in the thermodynamic limit. 

In the ferromagnetic case, close and below the critical point, kinetic energy 
fluctuations are physical and are due to the second-order phase transition. The LLE 
is related to these fluctuations and does not go to zero in the thermodynamic limit. 

In fig. 4 we analyze the behavior of the LLE in the limit of small energies for 
the antiferromagnetic case. As we have said previously the equations of the HMF 
are formally equivalent to those of a perturbed pendulum and in general we expect 
the HMF to be integrable when U tends to the ground state energy (i.e. when the 
perturbation goes to zero). This has already been checked for the ferromagnetic 
caseU" i and the law Ai oc U^^'^ , U ^0 was found; here we present the antiferro- 
magnetic case. We show in fig.4 Ai and Sks/^ vs. U* = U + 1/2 in log-log scale 
(the ground state energy in the repulsive case is —1/2). The dashed line indicates 
the presence of a power law with the same exponent 1/2 , i.e._Ai ~ (U*)^^"^ , U ^ 0. 
The same exponent was found also for nuclear-like systems Thus, also in this 
case there seems to be a universal law. However, though some heuristic arguments 
have been presented the deep theoretical reason of this law is not clear. 

§5. Conclusions 

We have discussed the most important results recently obtained for the HMF 
model. This model has revealed very useful for studying chaos in a Hamiltonian 
system with many degrees of freedom and in particular for undestanding the con- 
nection between microscopic chaos and macroscopic laws, i.e. thermodynamics. We 
have also discussed in particular the behavior of the Lyapunov exponents and the 
Kolmogorov-Sinai entropy for the ferromagnetic and antiferromagnetic case. While 
in the former case, where a second-order phase transition is present, one observes a 
well defined peak in the chaoticity indicators (Ai and Sks/^) a-t the critical point, 
in the latter case one has a smoother increase between the two integrable limits of 
very small and very large energy. In the high energy phase of both the ferromagnetic 
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Fig. 4. We show Ai and Sks/N as a function of U* — U + 0.5 for tiie antiferromagnetic case in 
log-log scale. Points are microcanonical numerical simulations for N=100. Dashed lines are fits 
of the behavior for very small energies. See text for more details. 

and the antiferromagnetic model Ai vanishes as N^^^^. In the ferromagnetic case 
the peak in Ai persists for — > oo and the LLE remains finite in the whole low 
temperature phase. Chaos persist in the thermodynamic limit. 

Though the border from low-dimensional to more realistic dynamical systems 
has been crossed and a lot of work has been done in this direction with the help 
of more powerful computers, we still have a long way to go in order to understand 
some important issues at the foundations of Statistical Mechanics. One has to admit 
that chaos in systems with many degrees of freedom is still poorly understood and 
represents the real challenge for the next decade. 
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